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ABSTRACT 



Context. X-ray observations provide a very good way to reveal the population of young stars in star forming regions avoiding the 
biases introduced when selecting samples based on infrared excess. 

Aims. The aim of this study was to find an explanation for the remarkable morphology of the central part of the S254-S258 star 
forming complex, where a dense embedded cluster of very young stellar objects (S255-IR) is sandwiched between the two HII 
regions S255 and S257. This interesting configuration had led to different speculations such as dynamical ejection of the B-stars 
from the central cluster or triggered star formation in a cloud that was swept up in the collision zone between the two expanding 
HII regions. The presence or absence, and the spatial distribution of low-mass stars associated with these B-stars can discriminate 
between the possible scenarios. 

Methods. We performed a deep Chandra X-ray observation of the S254-S258 region in order to efficiently discriminate young stars 
(with and without circumstellar matter) from the numerous older field stars in the area. 

Results. We detected 364 X-ray point sources in a 17' x 17' field (« 8 x 8 pc). This X-ray catalog provides, for the first time, 
a complete sample of all young stars in the region down to ~ 0.5 M B . A clustering analysis identifies three significant clusters: 
the central embedded cluster S255-IR and two smaller clusterings in S256 and S258. Sixty-four X-ray sources can be classified as 
members in one of these clusters. After accounting for X-ray background contaminants, this implies that about 250 X-ray sources 
constitute a widely scattered population of young stars, distributed over the full field-of-view of our X-ray image. This distributed 
young stellar population is considerably larger than the previously known number of non-clustered young stars selected by infrared 
excesses. Comparison of the X-ray luminosity function with that of the Orion Nebula Cluster suggests a total population of ~ 2000 
young stars in the observed part of the S254-S258 region. 

Conclusions. The observed number of ~ 250 X-ray detected distributed young stars agrees well with the expectation for the low- 
mass population associated to the B-stars in S255 and S257 as predicted by an IMF extrapolation. These results are consistent with the 
scenario that these two B-stars represent an earlier stellar population and that their expanding HII regions have swept up the central 
cloud and trigger star formation (i.e. the central embedded cluster S255-IR) therein. 

Key words. Stars: formation, low-mass, pre-main-sequence - X-ray: stars - Galaxy: clusters: individual: S254-S258 



1. Introduction 



,^ ; 1.1. The S254-S258 complex 



outflo ws dSnell & Ballvlll986tlMiralles et al.lll997fc Minier et al 
2005, l2007t iGoddi et al.l l2007b iJiang et alJ l200a IWang etal 



The south-eastern part of the molecular cloud complex in the 
Gem OBI association contains an embedded star-f orming re- 
gion with several diffuse HII regions (S254-S258. ISharplessi 
1959, see Fig. 1) . The most prominen t of these HII regions, 
S255 and S257 dChopinet et al.l Il974l) . are both powered by 
B0 stars, have diameters of ~ 4', and a projected separa- 
tion of ~ 6'. Sandwiched right between them is a dense 
dusty molecular cloud filament ( S255-IR or S255-2: lHever et all 
1989; D i Francesco et al.l 120081) that contains numerous em- 



bedded infrared sourc es ( Zinnecker etal J 1993 : Howard etal J 
19971: lltohet all 12001b lLongmore et alJl2006t lOiha et alJ l2006: 



Chavarria et al. 2008). Masers, HH-objects, jets, and molecular 



2011) provide clear evidence of very recent and ongoing star 
formation activity in t his cloud. The com bination of Spitzer mid- 
infrared observations (lAllen et al.ll2005l) with near-infrared im- 
ages of a 26' x 20' region l ed to the detection of 510 sources 
with near- or mid-IR excess (IChavarrfa et alj|2008l) . 87 and 165 
of which were classified as Class I and Class II sources, respec- 
tively. The large majority (80%) of these infrared excess-selected 
young stellar objects (YSOs) were found to be clustered. The 
central cluster S255-IR is the richest of these. It contains at leasQ 
140 infrared excess sources, among them 23 Class I sources. 
Another large fraction of the known YSO population is located 
in an elongated cluster at the southern edge of S256, and there 
are several smaller c lusters at different locations (see Fig. 12 in 
IChavarrfa et al . 200§. 



* Tables 2, 3, and 5 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



1 The census is incomplete because the Spitzer images of the dense 
cluster suffer from source crowding and saturation effects. 
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About 1' north of the center of S255 -IR, a strong far - 
infrared source, S255-N, was detected by iJaffe et all ([1984). 
This object is also associated with massive star formation 
and it is belie ved to be at an even earlier evolutionary stage 
than S255-IR dKurtz et alJ ll99l 12004: Cv eanowski et alJl2007t 



Wan g et all 201 lb . Another deeply embedded region, S255-S 



(Wan g et al. 



201 lb . is located about 1' south-west from S255- 



IR. It exhibits strong mm continuum emission (Wang et al. 201 1 ; 
iDi Francesco et ai1 l2008) but no other sign of active star fo rma- 
tion in near- and mid-IR observations. Min ier et"aD ([2007 ) sug- 
gested that this sub-region is in a very early pre-stellar phase 
of evolution. Scenarios for the spatial and temporal sequence 
of star formati on in the S254 - S258 c omplex have been rec ently 
discussed in IChayarria et al.1 d2008h . iBieging et all d2009b . and 
Wan g1it"aTl(l2011b . 

The distance of t he S254-S258 complex was only poorly 
known until recently. Pismis ~& Hassd (1 1976b and iMoffat et al.l 
(fl979h derived a v alue of 2.5 kpc; similar values (e.g., 
IChavarna et al.l 12008b . but also lower values down to 1 .5 kpc 
were used in later studies. iRvgl et al.1 (12010b performed high- 
precision astrometry using the 6.7 GHz methanol maser emis- 
sion from the source J0613+1708 in S255-IR and derived a very 
accurate trigonometric parallax of (1.59 + 0.07) kpc for S255. 
We will use this new and reliable distance for our present study. 

1.2. Evolution scenarios for the central region S255/257 

The most remarkable part of the S254-S258 complex and the fo- 
cus of the study presented here is the central region around the 
two HII regions S255 and S257. A very interesting, as yet un- 
explained, feature of this region is that the two B0 stars exciting 
these HII regions, ALS 19 and HD 253327, appear to be more 
or less iso lated and do not have any obvious co-spatial low-mass 
clusters (IZinnecker et al.lfl993b . This is rem arkable becaus e, ac- 
cording to the standard field star IMF (e.g., Kroupa 2001) each 
B0 star (M, ~ 15 M Q ) should be accompanied by about 300 
lower-mass stars. Fundamentally different possible explanations 
for the apparent absence of low-mass stars around these B0 stars 
have been proposed over the years. 

One scenario is based on the assumption of bimodal star 
formation. It assumes that the two high-mass B0 stars formed 
independently from the low-mass young stars in the central 
cluster and in a fundamentally different processes. The lack of 
clusters of low-mass stars around these two B0 stars would 
then imply that these high-mass stars formed in isolation (see 
Zinnecke r et alJI 1993b . One problem with this scenario is that in 
basically all other well investigated star forming regions, high- 
mass stars are always a s sociated with large n umbers of low-mass 
stars dTesti et al.lll999t iBricefio etai1l2007h . The case of S255 
and S257 would then represent a quite unique exception if the 
absence of low-mass stars was confirmed. 

The second scenario assumes that the two B0 stars formed 
together with the low-mass stars in the dense central clus- 
ter S255-IR, but were dynamically ejected, e.g. by means of 
close stellar encounters and N-body interactions. This model 
predicts that both B0 stars should move away from the cen- 
tral clusters with substantial velocities. Unfortunately, the avail- 
able HIPPARCOS proper-motions for the stars are not accurate 
enough to either support or rule out this scenario. 

A third scenario assumes multiple stellar generations and 
triggered star formation. Here, the two B0 stars belong to an 
earlier generation of stars that formed several Myr ago in this 
area. The expanding H II regions swept up diffuse gas and dust 
in their surroundings into shells, and formed the dense cloud in 



the interaction zone between them. This process of creating new 
clouds at the edges of shells or bubbles driven by high-mass stars 
is well established and observed at many locations (see, e.g . , 
iBrand et alJl201ll IZavagno etai]l20ld: iDeharveng etai] l2009). 
The particularly strong compression of the cloud at the intersec- 
tion of the two shells, caused by the ongoing expansion of the 
H II regions, may have triggered the formation of a new genera- 
tion of stars, i.e. the embedded cluster of young stellar objects. 
The ongoing expansion of the interacting bubbles would also 
provide a natural explanation why the youngest regions S255-N 
and S255-S are found just above and below the central young 
cluster S255-IR at the intersection of the shells. 

1.3. Importance of the low-mass stellar population 

A discriminant between the different evolution scenarios for 
the S255/257 region is the presence or absence of low-mass 
stars associated with the two B0 stars. While no low-mass stars 
should be present in the case of the bimodal star formation 
model or the dynamical ejection model, the multi-generation 
model predicts the presence of several hundred low-mass stars 
near these B0 stars, since the stellar populations in basically all 
well-i nvestigated OB ass ociations follow the standard field star 
IMF dBriceno et al. I l2007h . A handful emission line stars and a 
few dozen infrared excess objects (see lChavarria et al]|2008h are 
known inside or near the two H II regions, but their numbers 
are far to small for the expected low-mass population associated 
with the massive B0 stars. 

The apparent lack of associated low-mass stars may, how- 
ever, just be a result of the sensitivity limits of the existing ob- 
servations: a population of several Myr old low-mass stars would 
be quite hard to identify in the present optical and infrared im- 
ages, for several reasons. First, the low-mass stars would not be 
densely clustered around the B0 stars but scattered over a rather 
wide area, up to ~ 10 pc away from the massive stars, as typical 
for subgroups in OB associations. Second, these low-mass stars 
would be quite hard to see in most existing optical or infrared 
images of the region: since a ~ 10 Myr old 1M [0.2M o ] stars 
should have magnitudes of V > 18.3 [22.2] and K > 14.7 [17.0], 
they would be not easy to detect in the nebulosity of the H II re- 
gions and the diffuse infrared emission in this region. Third, even 
the availability of Spitzer observations is of limited use here: al- 
though Spitzer data are sensitive enough to detect a good frac- 
tion of the low-mass stars, the Spitzer images of the region are 
dominated by unrelated field stars (note that this region lies very 
close to the galactic plane). The usual approach to identify young 
stars by their infrared excesses is not feasible here, because at 
an age of more than a few Myr, most of the low-mass associ- 
ation members have already lost their cir cumstellar disks an d 
thus should not exhibit infrared excesses dBriceno et al. I l2007h . 
It is thus impossible to identify and distinguish a population of 
several (~ 3 — 10) Myr old low-mass association members from 
unrelated field stars with optical or infrared photometry alone. 

Sensitive X-ray observations can provide a very good solu- 
tion of this problem, since they allo w to detect the young stars 
by their strong X-ray emission (e.g., Feige lson et alJ[20071) and 
efficiently discriminate them from the numerous older field stars 
in the survey area. The median X-ray luminosity of % 10 Myr 
old solar-mass stars is ss 10 ergs -1 ; th is is nearly 1000 times 
highe r than for solar-mass field stars (see lPreibisch & Feigelsonl 
2005), and makes these young stars relatively easily detectable 
X-ray sources. Another very important aspect is that X-ray ob- 
servations trace magnetic activity rather than photospheric or cir- 
cumstellar disk emission from young stars, and are thus comple- 
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mentary to the available optical and infrared data of the region. 
The X-ray selected sample of low-mass stars will be not biased 
toward stars with circumstellar disks identified in the Spitzer 
data. Furthermore, an X-ray image is not subject to confusion 
from bright diffuse emission by heated gas and dust. X-rays 
can penetrate deeply into obscurin g material and are ver y effec- 
tive in detecting embedded YSOs ( Getma n et al ] l2005ah . Many 
X-ray studies of star forming re gions have demonstrated the 
success of this method (see, e.g.. | Preibisch & Zinnec ker 2002; 
Broos et all l2007t iForbrich & Preibischl 120071: iTownslev et al. 



201 lb . Also, the relations between the X-ray properties and ba- 



sic stellar properties in young stellar populations are now very 
well established from very deep X-ray observatio ns such as the 
Chandra Orion Ultradeep P roject (COUP) (see [Getman et al. 
2005b; PreibischetaLl|2005). To summarize, a deep X-ray im- 
age of the S254-S258 complex can reveal the full young stellar 
populations in the area and provide essential information about 
the star formation history. 

At distances beyond 1 kpc, very good angular resolution is 
required to resolve the individual sources in the dense young 
clusters and to allow a reliable identification of the X-ray sources 
with the numerous infrared sources (note that the complex is 
almost exactly on the galactic plane, b = -0.048°). The Chandra 
X-ray observatory, that provides an on-axis PSF of < 1", is the 
only currently active X-ray mission that has sufficient angular 
resolution for this purpose. 

We have therefore performed a deep Chandra X-ray ob- 
servation of this extraordinary star forming region in order to 
uncover the population of low-mass association members. Our 
study focuses on the central region of the S254-S258 complex, 
i.e. the two H II regions S255 and S257 and the embedded cluster 
S255-IR between them. A characterization of the size, the spatial 
distribution, and the properties of the low-mass population can 
provide important information on the star formation history and 
discern between the different models for the relation between 
S255/S257 and the embedded cluster S255-IR. In Section |2] we 
describe the Chandra observations and data reduction. Section 
3 presents the basic X-ray properties of the detected sources. 
Section |4] analyzes the X-ray population of the S254-S258 star 
forming complex, and Section 5 discusses the spatial distribution 
of the X-ray sources and the implications on the star formation 
process in S254-S258. A more detailed analysis of the optical 
and infrared properties of the individual X-ray detected young 
stars (that can provide direct information on the ages, masses, 
and the circumstellar disks around these stars) will be presented 
in a forthcoming paper. 



2. Observations and data reduction 

2.1. Chandra X-ray observations of S254-S258 

The S254-S258 complex was observed (PI: Th. Preibisch) 
in November 2009 with the Imaging Array of the Chandra 
Advanced CCD Imaging Spectrometer (ACIS-I). ACIS-I pro- 
vides a field of view of 17' x 17' on the sky. At the 1.6 kpc 
distance of S254-S258 this corresponds to 7.9 x 7.9 pc. The 
aimpoint of the observation was 0^2000) = 06 h 12 m 54.0 s , 
6 = +17° 59' 24". The observation was performed in the stan- 
dard "Timed Event, Faint" mode (with 3x3 pixel event islands). 
The total net exposure time of 74725 s (20.76 hours) was split 
into two parts, separated by about 4 days. The details of these 
two observation parts are given in Table Q] Since the roll an- 
gles (i.e. the orientation of the detector on the sky) is equal for 
both observations, each source is at the same detector position in 



both parts, making merging of these two data sets rather straight- 
forward. Two of the CCDs of the ACIS-S spectroscopic array 
were also turned on during our observations. However, since the 
PSF at the large off-axis angles at these detectors is strongly de- 
graded, their point-source sensitivity is reduced; only six X-ray 
sources are detected in the field of these ACIS-S chips. 

The basic data products of our observation are the two Level 
2 processed event list provided by the pipeline processing at 
the Chandra X-ray Center, that list the arrival time, location 
on the detector and energy for each of the 626 140 detected X- 
ray photons. We combined the two pointings with the merge jail 
script, a Chandra contributed software that make use of standard 
tools. The mean background count rate in our merged 
image, determined from several large source-free regions, is 
2.95 x 10~ 7 counts s _1 pixeF 1 , corresponding to a mean back- 
ground level of 0.02 counts pixel 1 . 

At a distance of 1.6 kpc, the expected ACIS point source 
sensitivity limit for a 5-count detection on-axis in a 75 ks ob- 
servation is Lx.min ~ 10 295 erg s _1 , assuming an extinction of 
Ay < 2.5 mag (Nh < 5 x 10 21 cirT 2 ) as typical for the stars 
in the HII regions, and a thermal plasma wit h kT — 1 keV 
(whic h is a typical value for young stars; see, e.g JPreibisch et ail 
120051) . Using the empirical relation between X-ray luminosity 
and stellar mass and the temporal evolution of X-ray luminosity 
from the sample of young stars in the Orion Nebula Cluster that 
was very well studied in the Chandra Orion Ultra deep Project 
dPreibisch et al.ll2005t IPreibisch & Feigelsonll2005l) . we can ex- 
pect to detect almost all stars in S254-S258 with masses greater 
than 0.5 M and about half of the 0.1 - 0.5 M stars. The ex- 
pected level of detection completeness is ^ 90% for stars with 
M, > 0.5 M Q (corresponding approximately to spectral types 
earlier than Ml) and drops below 50% at M„ < 0.25 M (spec- 
tral types SS M5). Note that these values are valid for the central 
part of the observed field; sensitivity is ~ 3 - 4 times worse at 
the edges of the ACIS field. 

2.2. Source detection and X-ray source catalog 

The source detection was performed in a two-step process. The 
first detection step was performed in a rather aggressive man- 
ner in order to find even the weakest possible sources, deliber- 
ately accepting some degree of false detections. In the second 
step, this list of potential sources was then cleaned from spuri- 
ous detections by a deta iled individual analy sis. We employed 
the wavdetect algorithm (Fre eman et al.ll2002l a CIAO mexican- 
hat wavelet source detection tool) for locating X-ray sources in 
our merged image, and used a rather low detection threshold of 
10~ 5 . This step was performed in three different energy bands, 
the total band [0.5-8.0] keV, the soft band [0.5-2.0] keV, and the 
hard band [2.0-8.0] keV, and with wavelet scales between 1 and 
16 pixels. We also performed a visual inspection of the images 
and added some 30 additional candidates to the merged catalog 
from the wavelet analysis, resulting in a final catalog of 5 1 1 po- 
tential X-ray sources. 

To clean this catalog from spurious sources, we then 
performed a detailed analysis of each individual candidate 
source with the AC IS Extract (AE hereafter) software packag^l 
dBroos et al.l l2010h . A f ull description of the p rocedures used 
in AE can be found in [Getman et al.l d2005bl) . Townslev et al.l 
(2003]) and iBroos et al.l d2007l) . The following three steps were 



2 Chandra Interactive Analysis of Observations, version 4.2: 
http://cxc.harvard.edu/ciao/index.html 

3 http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 
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Fig. 1. Top: Negative grayscale representation 
of the optical image of the S254-S258 complex 
from the Digitized Sky Survey. The black el- 
lipses represent the five H II regions that define 
the complex. The location of the central, op- 
tically invisible embedded cluster S255-IR is 
marked by the arrow. 

Center: Spitzer IRAC 4 image of the cen- 
tral part of the S254-S258 complex. This im- 
age was created from the basic calibrated data 
products for the programs 201 and 30784 re- 
trieved from the Spitzer archive and mosaicked 
with the MOPEX software available from the 
Spitzer Science Center. Note that parts of the 
bright emission from the central embedded 
cluster S255-IR is saturated in these data. 
Bottom: Chandra ACIS-I image of S254-S258 
in the [0.5-8.0] keV band. Blue ellipsoids rep- 
resent extraction regions for the individual de- 
tected X-ray sources based on a model of the 
local PSF that encircles 90% of total energy. 
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Obs.Id. 


Date 


Start - End Time [UT] 


Exposure time 


Level 2 events 


10983 


2009-11-16 


11:40:48-23:35:08 


40570 s 


340411 


12022 


2009-11-20 


05:14:15 - 15:18:20 


34 155 s 


285 729 



performed by AE in order to prune our input catalog from spuri- 
ous detection (including afterglows): 

1 . Extraction regions were defined as the 90% contours of the 
local PSF (or smaller in the case of other nearby sources), 
and source events were extracted. Energy dependent correc- 
tions for the finite extraction regions were applied; 

2. Local background events were extracted after masking all 
the sources in the catalog; 

3. The Poisson probability (Pb) associated with the "null hy- 
pothesis", i.e. that no source exist and the extracted events 
are solely due to Poisson fluctuations in the local back- 
ground, is computed for each source. 

All candidate sources with Pb > 0.01 were rejected as back- 
ground fluctuations. After 8 iterations of this pruning procedure 
our final catalog consisted of 364 sources. It contains 344 pri- 
mary sources with Pb < 0.003, and 20 tentative sources with 
0.003 < Pb < 0.01. The extraction regions for the sources in our 
final catalog are plotted on the Chandra image in Fig.Q] 

2.3. X-ray point-source analysis with ACIS Extract 

The AE software also determines basic properties for each of the 
detected sources, such as the net (i.e., background-subtracted) 
counts in various energy bands, the median photon energy, sta- 
tistical test for variability, and a measure of the incident pho- 
ton flux. These properties are reported in Table 2 (available in 
the electronic edition). Sources are sorted by increasing right 
ascension and identified by their sequence number (Col. 1) or 
their IAU designation (Col. 2). While the general X-ray proper- 
ties were determined from the merged data set (the AE software 
is well suited for this purpose), we note that the spectra (see 
Sect. 13.1.31 were extracted from the individual observations. 

2.4. Expected contamination of the X-ray source sample 

As in any X-ray observation, there must be some degree of 
contamination by galactic field stars as well as extragalactic 
sources. To quantify the expected level of this contamination, we 
consider the results from the recent Chan dra Carina Complex 
Project (CCCP; see iTownslev et all 1201 ll) . for which the indi- 
vidual pointings had very similar exposure times (* 60 - 80 ks) 
as our S254-S258 pointing. Furthermore, S254-S258 is at nearly 
the same galactic latitude as the Carina Nebula, suggesting that 
the background contamination should be very similar in these 
two regions. 

Fo r the CCCP data set, the classification study of lBroos et alj 
(12011 ah . which considered the X-ray, optical, and infrared prop- 
erties of the sources (that differ for the different contaminant 
classes), found that 716 X-ray sources in the 1.46 square-degree 
CCCP survey are are foreground stars, 16 are background stars, 
and 877 are extragalactic (AGN) contaminants. Scaling these 
numbers to the field-of-view of our S255 pointing gives 39 fore- 
ground stars, 1 background star, and 48 extragalactic (AGN) 
contaminants. However, since S254-S258 is considerably closer 
(1.6 kpc) than the Carina Nebula (2.3 kpc), the number of fore- 
ground stars should be accordingly smaller, approximately by a 



factor of (1.6/2.3) 3 ~ 0.34. Furthermore, the number of fore- 
ground stars in the Carina Nebula is particularly high since this 
direction is close to the tangent point of the Carina spiral arm. 
These considerations imply that the contamination in our S254- 
S258 field should be clearly dominated by ~ 48 expected extra- 
galactic sources (AGNs). 

A characteristic of extragalactic contaminants is that their 
optical and infrared counterparts should be very faint, in fact 
mostly undetected in the available optical and infrared images. 
As we describe in more detail below, our search for counterparts 
of the X-ray sources left 46 X-ray sources outside the central 
embedded cluster S255-IR without optical or infrared counter- 
parts. This number agrees well with the expected number of ex- 
tragalactic contaminants. 

Assuming at most 10 contaminating foreground/background 
stars, the total expected number of contaminants would be S5 58. 
With a total number of 364 detected X-ray sources, the expected 
level of contamination for our sample is thus S= 15%. 

3. Properties of the X-ray source in S254-S258 

3.1. X-ray fluxes and luminosities 

An accurate determination of the intrinsic X-ray source lu- 
minosities requires good knowledge of the X-ray spectrum. 
However, for the majority of the X-ray sources the number of 
detected photons is too low for a detailed spectral analysis. Only 
25 sources in our catalog have more than 80 net counts, the 
practical lower limit for meaningful spectral analysis. For these 
bright sources we performed a detailed spectral fitting analysis 
to derive the plasma temperature and the extinction, and from 
these quantities we can calculate the intrinsic (i.e. , extinction- 
corrected) X-ray luminosities, as described in detail in Section 

E3 

For the weaker X-ray sources, for which a meaningful spec- 
tral analysis is not feasible, one cannot determine intrinsic X-ray 
luminosities without knowledge of the extinction. This is a sub- 
stantial problem because the young stars in the S254-S258 com- 
plex show a very wide range of extinctions. There are numerous 
optically visible stars with low obscuration (at most a few mag- 
nitudes of visual extinction), while other stars suffer from cloud 
extinction up to about Ay ~ 20 mag, and embedded YSOs show 
additi onal circumstellar ext inctions up to Ay ~ 50 mag and be- 
yond (Chavam'a et al. 2008). This implies that we cannot simply 
use a common count-rate to flux conversion factor to determine 
intrinsic X-ray luminosities but have to consider each source in- 
dividually. 

3.1 .1 . Observed X-ray fluxes 

An estimate of the observed (i.e. not the intrinsic) X-ray flux 
is computed by AE. This quantity, called FLUX2, is calculated 
from the number of detected photons and using a mean value of 
the instrumental effective area (through the Ancillary Response 
Function, ARF) over energy. The FLUX2 values (derived for the 
full band, i.e. [0.5 - 8] keV range), are reported in column (3) in 
Table 3 (available in the electronic edition). It should be noted 
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Fig. 2. The blue histogram shows the distribution of FLUX2 val- 
ues for the X-ray sources in S254-S258. For co mparison, the 
distribution of fluxes from the COUP data from Getman et al. 
d2005bl) is shown by the black histogram, scaled to the distance 
of S254-S258, i.e. 1.6 kpc. 



3.1.2. X-ray luminosities from XPHOT 

An estimate of the intrinsic, i.e. extinction corrected, X-ray lu- 
minosity for sources that are too weak for a detailed spectral 
an alysis can be obtained with the XPHOT softwarfl developed 
bv iGetman et alJ (l2Q 1 Oh . XPHOT is based on a non-parametric 
method for the calculation of fluxes and absorbing X-ray col- 
umn densities of weak X-ray sources. X-ray extinction and in- 
trinsic flux are estimated from the comparison of the apparent 
median energy of the source photons and apparent source flux 
with those of high signal-to-noise spectra that were simulated 
using spectral models characteristic of much brighter sources of 
similar class previously studied in detail. This method requires 
at least 4 net counts per source (in order to determine a meaning- 
ful value for the median energy) and can thus be applied to 255 
of our 364 sources. Columns (4) to (7) of Table 2 report appar- 
ent and intrinsic (corrected for absorption, noted with subscript 
c) luminosities in the hard and total band, assuming a distance 
of 1.6 kpc. The resulting intrinsic X-ray luminosities range from 

1Q 29.4 tQ 10 32.3 erg s -l 

Figure [3] shows the distribution of median photon ener- 
gies and the deduced hydrogen column densities estimated by 
XPHOT. The median value of the derived hydrogen colum den- 
sities is log^H [cm -2 ]) = 22.04, corresponding to a visual ab- 
sorption of Ay ~ 6 magnitudes. If we consider the sub-sample 
of sources located in the central embedded cluster S255-IR, this 
value rises to log (/V# [ctrT 2 ]) = 22.21 (Ay ~ 9 magnitudes), 
clearly showing stronger obscuration for the embedded sources. 



that this FLUX2 values suffer from a systematic error with re- 
spect to the true incident flux, because the use of a mean ARF 
is only correct in the hypothetical case of a flat incident spec- 
trum, an assumption that probably not fulfilled. Nevertheless, the 
FLUX2 represents the best flux estimate that can be obtained for 
weak sources. 

In Figure [2] we compare the distribution of FLUX2 values 
for the S254-S258 sample to that of the stars in the Orion Nebula 
Cluster obtained in the context of the Cha ndra Orion Ultradeep 
Project (COUP, see iGetman et alJl2005bl) . Note that the COUP 
fluxes were scaled to the 1.6 kpc distance of S254-S258. The 
two distribution show a very similar shape in the range between 
» 10~ 60 photons cirT 2 s _1 and m 10~ 4 5 photons ctrT 2 s _1 . The 
differences between the two distributions can be explained as 
follows: 

First, the COUP sample shows a few stars with fluxes of 
> 10~ 45 photons cirT 2 s , while no such very bright sources 
are seen in S254-S258. These very bright sources are the highly 
X-ray luminous O-type stars in the Orion Nebula Cluster. The 
fact that the S254-S258 does not contain such high-mass stars 
explains the absence of similarly high values in the observed 
distribution of incident fluxes for S254-S258. 

Second, the peak and turn-over of the S254-S258 distribu- 
tion at as 10~ 6 ' 2 photons cirT 2 s is a direct consequence of 
the higher sensitivity limit of our S254-S258 X-ray observation. 
As S254-S258 is about 4 times more distant than the ONC, and 
since the exposure time of our S254-S258 Chandra observation 
is less than one tenth of the 840 ks COUP observation, the ex- 
pected sensitivity limit should be about 150 times higher. 

Third, the number of sources per bin is always lower for 
S254-S258 compared to the ONC. This suggests that the total 
stellar population in the observed part of S254-S258 is smaller 
than in the ONC as observed in the COUP. 



3.1 .3. X-ray spectral fits of bright sources 

For the 25 sources in our sample with more than 80 net counts 
we performed a s pectral fi t ting a nalysis using AE and the XSPEC 
software vl2.5 dArnaudl 1 19961) . We used model s with one- 
or tw o-temperature thermal VAPEC components ( Smit heTaTI 
1200 ll) and the TBABS multiplicative model to describe the ef- 
fect of extinction by interstellar (and circumstellar) material (as 
measured by the hydrogen column density Nh)- The plasma 
abundances for the VAPEC comp onents were fixed at the val- 
ues adopted by the XEST study dGttdel et alj|2007l) to be typi- 
cal for pre-main sequence staro For the extinction we used the 
sta ndard interst ellar abundances in the TBABS model as listed 
in IWilms et alJ d2000l) . In order to evaluate the goodness of 
our fits we c hoose t o app l y the C-statistic (a m aximum likeli- 
hood method; ICashl (fl979T) : IWachter et al.l (0971 1 which is bet- 
ter suited than the classic x 1 statistic for low-count data. 

For two sources, # 97 and # 230, a two-temperature model 
was required for an acceptable fit. The remaining spectra are well 
fit with a single thermal component. A few selected examples of 
the spectral fits are shown in Figure [4] The spectral parameters 
are reported in Table [4] The hydrogen column densities range 
from log (jVtf [cirT 2 ]) = 20 to log (n h [cirT 2 ]) = 23.08, cor- 
responding to visual absorptions between Ay ^ 0.1 mag and 
Ay ~ 65 mag. These values are in good agreement with the es- 
timates derived with XPHOT. The median value is 22.04, corre- 
sponding to Ay ~ 6 mag. The plasma temperatures range from 
~ 0.5 keV (6 MK) up to ~ 15 keV (170 MK). In Tableg]we also 
report luminosities derived from the spectral fit, assuming a dis- 

4 http://www.astro.psu.edu/users/gkosta/XPHOT/ 

5 The adopted abundances, relative to the solar photospheric abun- 
dances given by And ers & Grevessd {1989), are: C = 0.45, N = 0.788, 
O = 0.426, Ne = 0.832, Mg = 0.263, Al = 0.5, Si = 0.309, S = 0.417, 
Ar = 0.55, Ca = 0.195, Fe = 0.195, Ni = 0.195. 
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Fig. 3. Distributions of median photon energies (left) and hydrogen column densities (right) determined by XPHOT (right) for the 
whole sample (black histograms). The blue dashed histograms show the distributions restricted to the sub-sample of sources located 
in the 1 arcmin radius region centered on the embedded cluster S255-IR. 



tance of 1.6 kpc. The intrinsic luminosities are calculated from 
the spectral fit parameters, setting extinction to zero. The range 
of extinction-corrected intrinsic X-ray luminosities spans from 
logL v = 30.56 ergs" 1 to 32.11 erg s" 1 for the full [0.5-8.0] keV 
band. 

The X-ray properties can give us some clues about the nature 
of the sources. The majority of the X-ray sources has plasma 
temperatures and X-ray luminosities in the typic al ranges found 
for Y SOs in other star forming regions (see, e.g., Prei bisch et alJ 
120051) : together with the fact that most X-ray sources have 
clear optical/infrared counterparts, this suggests that these X- 
ray sources actually are young stars in the S254-S258 complex. 
However, the sources for which the spectral fit yields extremely 
high plasma temperatures of kT > 7.5 keV deserve special at- 
tention. Such very hard spectra are typical for extragalactic ob- 
jects (AGNs), but also someti mes found for very y oung stellar 
objects (protostars) (see, e.g.. Ilmanishi et al] |2001h . As proto- 
stars are usually deeply embedded in the clouds in which they 
formed, the location of these very hard X-ray sources provides 
another hint towards their likely nature. 

Considering these issues, source number 2 may well be an 
extragalactic contaminant. First, the spectral fit yields an ex- 
tremely high temperature, but only moderate extinction. Second, 
this source has no optical/infrared counterpart, and, third, it is 
located at the periphery of the S254-S258 region, well outside 
the boundaries of the molecular clouds. Furthermore, its X-ray 
spectrum can also be well reproduced with a power-law model, 
as typical for AGN X-ray sources. Similar arguments apply to 
source 30. 

The other X-ray sources with extremely high plasma tem- 
peratures are located in dense clouds, have infrared counterparts, 
and are thus probably deeply embedded very young stellar ob- 
jects. 

It is interesting to compare the X-ray luminosities from the 
spectral fits to those derived with XPHOT. Figure [5] shows that 
the results from these two different methods agree quite well 
for the majority of cases. Only for four of the highest luminos- 
ity objects we find that XPHOT seems to systematically over- 




31 32 
>°g V-phot [erg/s] 



33 



Fig. 5. Comparison of the intrinsic full band [0.5-8.0] keV lu- 
minosities derived from the spectral fits to the 25 bright sources 
(TableHJl to the intrinsic full band [0.5-8.0] keV luminosities de- 
termined with XPHOT. 



estimate the X-ray luminosities by * 0.1 - 0.2 dex. The gener- 
ally good agreement suggests that the X-ray luminosities derived 
with XPHOT are reliable. We therefore will use the XPHOT re- 
sults for our analysis of the X-ray luminosity function presented 
below. 
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Fig. 4. Chandra X-ray spectra and best-fit models of four bright X-ray sources in S254-258. The crosses show the measured spectra, 
the solid lines show the best-fit models; in the two cases where a two temperature model was required, the dotted lines show the two 
individual spectral model components. CXOU J061244. 17+175914.0 is the B0.5 star HD 253327 that illuminates the S257 nebula. 
CXOU J061253. 73+175724.7 is the source with the highest derived X-ray luminosity in our sample; it is located in the embedded 
cluster S255-IR. CXOU 061255.04+175930.7 is another embedded infrared source in S255-IR. CXOU 061312.58+180706.2 is the 
source with the second highest X-ray luminosity in our sample. 



3.2. X-ray source variability 

A first analysis of the time variability of individual X-ray sources 
is performed by AE by comparing the arrival times of the indi- 
vidual source photons in each extraction region to a model as- 
suming temporal uniform count rates. The statistical significance 
for variability is quantified computing the 1 -sided Kolmogorov- 
Smirnov statistic (Col. 15 of Table 2). In our sample, 23 sources 
show significant X-ray variability (probability of being con- 
stant P const < 0.005) and 19 are classified as possibly variable 
(0.005 < 

'const 

< 0.05). 

The light curves of the variable X-ray sources show a variety 
of temporal behavior; six of the most interesting lightcurves are 
shown in Fig. [6] Five of these sources show flare-like variability, 
i.e. a fast increase of the count rate followed by a slow exponen- 
tial decay, as typical for solar-like magnetic reconnection flares 
(see, e.g., iFavata etalj[2005t IWolk et alj l200l iPreibisch et all 
fl993h . The other variable sources show irregular variability or 
slowly increasing or decreasing countrates. This kind of X- 
ray variability is t ypical for young stellar objects (see, e.g., 
IStassun et ai]|2006h . 

We also investigated time variability by comparing the count 
rate for each source of the two single Chandra pointings, i.e. at 



a time difference of about 4 days (see Figure |7). One can see 
that the majority of sources show changes in the count rates by 
less than a factor of 2. Only for 21 sources the count rates dif- 
fer by more than 3cr between the two observations. The result 
is consistent with the assumption that the X-ray emission from 
young stellar objects is a superposition of many flares of differ- 
ent amplitude, where weak flares are very frequent while very 
strong flares oc cur more rarely, at r ates of about one such event 
per week (e.g., Get man et alj |2008). 

4. Characteristics of the X-ray stellar population of 
S254-S258 

4. 1 . Optical and infrared counterparts of the X-ray sources 

In order to identify counterparts of the X-ray sources in other 
wavelengths, we used the optical images from the Digitized 
Sky Survey (DSS), the Two Micron All Sky Survey (2MASS) 
point source catalog , and the Spitzer-IRAC catalog from 
IChavarrfa et al.l (I2008I) . The results of the cross-correlation are 
reported in Table 5 (available in the electronic edition). Our 
visual inspection of the red and blue DSS images gave opti- 
cal counterparts to 95 X-ray sources (i.e. 26% of all 364 X-ray 
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Fig. 6. Lightcurves for six significantly variable sources. The solid dots show the arrival time (measured from the start of the 
observation) and the energy of each of the detected source photons. The histograms show the corresponding binned lightcurves. 
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Table 4. Spectral parameters of brighter sources: the spectral fit was performed with an absorbed thermal plasma model with one 
(tbabs*vapec) or two components (tbabs*(vapec+vapec)). 



Source 


CP/ 


CXOU J 


Net Counts 


log N H 


kT, 


kT 2 


log L s 


log Ui 


log L h c 


logL, 


log L, f 


No. 


IR class 




(counts) 


(cm"") 


(keV) 


(keV) 


(ergs" 1 ) 


(ergs -1 ) 


(erg s" 1 ) 


(ergs -1 ) 


(ergs Jl ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


ves 


061147.81 + 180312.6 


443.3 


21.81+0.15 


0.7+0. 1 




30.89 


30.14 


30.19 


30.97 


31.47 


2 


no 


061154.83+180016.5 


192.1 


22.04+0.46 


15.0+5.1 




30.26 


31.08 


31.11 


31.14 


31.27 


3 


yes 


061159.55+175802.9 


159.6 


22.38+0.64 


2.6+1.2 




30.21 


30.96 


31.05 


31.03 


31.40 


5 


yes / III 


061206.65+180336.6 


239.0 


21.41+0.16 


0.5+0.1 




30.88 


29.55 


29.57 


30.90 


31.17 


27 


ves 


061230.43+175506.8 


81.8 


20.48+0.19 


9.4+9.9 




30.23 


30.53 


30.53 


30.70 


30.71 


30 


ves 


061230.98+180336.9 


386.2 


22.07+0.13 


15.0* 




30.62 


31.46 


31.49 


31.52 


31.65 


32 


ves 


061231.17+180853.8 F 


99.5 


21.96+0.41 


3.8+2.3 




30.28 


30.78 


30.81 


30.90 


31.09 


97 


yes / III 


061244.17+175914.0 


174.6 


21.56+0.26 


0.6+0.2 


3.4+1.6 


30.57 


30.45 


30.46 


30.81 


31.02 


102 


ves 


061244.76+175946.8 


156.6 


21.61+0.15 


3.1+0.8 




30.41 


30.65 


30.67 


30.85 


30.98 


104 


ves 


061245.23+175810.3 F 


116.6 


20.0* 


1.5+0.2 




30.45 


30.06 


30.06 


30.60 


30.61 


113 


yes / III 


061246.71 + 175418.2 


86.6 


21.53+0.23 


0.7+0.1 




30.38 


29.47 


29.49 


30.43 


30.73 


147 


yes / III 


061251.04+180644.1 


141.8 


20.00+16.03 


1.0+0.1 




30.69 


29.66 


29.66 


30.73 


30.74 


183 


ves 


061253.48+175633.9 


91.2 


23.08+3.43 


54.2+724 




28.43 


31.13 


31.35 


31.13 


31.48 


190 


yes 


061253.73+175724.7 


373.2 


22.99+1.25 


5.4+1.2 


_ 


29.51 


31.64 


31.88 


31.65 


32.11 


209 


yes 


061254.33+175927.7 F 


86.5 


22.86+1.99 


8.6+10.8 




29.13 


30.99 


31.16 


31.00 


31.35 


230 


yes 


061255.04+175930.7 F 


250.5 


22.30+0.53 


0.9+0.4 


8.7+13.0 


30.44 


31.11 


31.18 


31.20 


31.61 


255 


yes / 1 


061258.21 + 175848.1 


91.7 


22.34+0.53 


2.6+0.9 




30.00 


30.70 


30.79 


30.78 


31.14 


333 


yes 


061312.58+180706.2 


379.7 


22.56+0.46 


7.9+3.5 




30.33 


31.60 


31.70 


31.63 


31.90 


337 


no 


061316.04+175416.7 


139.9 


22.22+0.42 


3.9+1.6 




30.19 


30.91 


30.96 


30.98 


31.24 


341 


yes 


061317.11 + 175344.9 F 


265.6 


21.74+0.17 


5.0+1.6 




30.70 


31.16 


31.17 


31.29 


31.41 


348 


yes / II 


061325.57+175230.5 


163.1 


22.18+0.39 


1.6+0.3 




30.44 


30.73 


30.81 


30.91 


31.35 


349 


yes / III 


061325.94+175923.5 F 


95.6 


20.0* 


1.0+0.3 




30.48 


29.73 


29.73 


30.55 


30.56 


354 


yes 


061327.41 + 175554.3 


125.2 


21.91+0.3 


5.6+4.1 




30.25 


30.82 


30.85 


30.93 


31.08 


355 


yes 


061327.61 + 175517.8 


182.8 


22.11+0.35 


3.8+1.4 




30.38 


31.00 


31.05 


31.10 


31.32 


359 


yes 


061328.37+175604.4 


97.7 


22.15+0.52 


4.3+2.6 




30.07 


30.75 


30.80 


30.84 


31.06 



* indicates frozen parameters in the fit. 

Col. (2): presence of an optical or infrared counterpart, infrared class from Spitzer photometry (if available). 

Col. (3): sources flagged with "F" showed flare-like variability during our observations; lightcurves are shown in Fig. [6] 

Col. (4): absorbing hydrogen column density of the best-fit. 

Cols. (5) and (6): plasma temperature(s) of the best-fit. 

Cols. (7) to (1 1): X-ray luminosities (for an assumed distance of 1.6 kpc) in the soft (s, [0.5 - 2.0] keV) band, the hard (h [2.0 - 8.0] keV) band, and the total (f [0.5 - 8.0] keV) band. 
Absorption-corrected luminosities are denoted with the subscript c. 



sources). Our cross-correlation with the 2MASS Point Source 
Catalog lead to 231 near-infrared counterparts (i.e. a counterpart 
fracti on of 63%). Our cross- correlation with the infrared catalog 
from Chavam a et al.l ((2008) yielded 293 infrared counterparts, 
i.e. 80% counterpart fraction. For 58 X-ray sources we did not 
find a counterpart in any of the inspected optical and infrared im- 
ages. Twelve of these sources are located in or very close to the 
central embedded cluster S255-IR. These X-ray sources may be 
very deeply embedded protostars or young stellar objects located 
behind the dense molecular cloud clumps; the non-detection of 
optical/infrared counterparts would then be related to very high 
extinction. The remaining 46 X-ray sources without known opti- 
cal/infrared counterpart outside this cluster show a rather homo- 
geneous spatial distribution, as expected for (mostly extragalac- 
tic) contaminants. 

It is interesting to consider the infrared classification of these 
sources based on th e IRAC spectral energy distribution (SED) 
slope determined by IChavarrfa et al. (2008). Unfortunately, the 
matching of our X-ray source-list with this infrared catalog is not 
straightforward. The catalog contains 26 821 infrared sources. 
However, most of these are only detected in the deep near- 
infrared images, and just about 6400 of these are detected in the 
Spitzer data. Infrared classifications are only available for the 
462 infrared sources that are detected in all four IRAC bands. 
The majority of the catalog entries are very faint NIR sources, 
and many of these are probably background objects rather than 
young stars in S254-S258. Furthermore, the fact that only 4225 
of the 26 821 sources are detected in all three of the J-, H-, and 
/T-bands also suggest that there may well be a significant number 
of spurious detections among the faint sources detected in only 
one band. This very high number of faint infrared sources pro- 
duces serious problems in any attempt to find the correct infrared 
counterparts for our X-ray sources, since many Chandra sources 



have more than one possible counterparts within the X-ray error 
radius. In these cases, the closest positional match is not nec- 
essarily the true counterpart. Due to the increasing number of 
infrared sources at fainter magnitudes, good positional matches 
with very faint infrared sources may in fact be just chance super- 
positions of physically unrelated sources, and one of the other 
possible matches may be the true counterpart A reliable iden- 
tification of the infrared counterparts requires a sophisticated 
approach and will be addressed in the next step of our study. 
Nevertheless, we can mention here the results of a preliminary 
source matching, where we only considered the spatially clos- 
est match to each X-ray source. We find that 8 X-ray sources 
have closest matches classified as Class I YSOs (embedded very 
young stellar objects with infalling envelopes), 50 X-ray sources 
have closest matches classified as Class II YSOs (Classical T- 
Tauri stars, CTTs), and 8 X-ray sources have closest matches 
classified as Class III YSOs ("Weak line" T-Tauri stars, WTTs) 
in the infrared catalog. 

4.2. The X-ray luminosity function 

The X-ray luminosity function (XLF) is the product of the 
distribution of X-ray luminosities of stars with a given mass 
and the number of stars per mass interval, i.e. the initial 
mass function (IMF). Although the correlation between stel- 
lar mass and X-ray luminosit y shows a considerable scatter 
(see, e.g., iPreibisch et al.ll2005l) . X-ray studies of a large num- 
ber of young stellar clusters have shown that the XLF appears 
to be rather universal and constant in different environments 



6 We note that similar problems were encountered in an X-ra y and 
infrared study of the Carina Nebula; see Preibisch et al. (2011) for a 
more detailed discussion. 
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Fig. 7. This plot shows for each X-ray source the net counts in 
the first pointing compared to the net counts in the second point- 
ing (obtained about 4 days later). The sequence of blue error 
bars at the upper and right edge show the Poisson statistical un- 
certainties for different numbers of net counts. The solid line 
indicates the expected relation for sources with constant count 
rates in the two pointings; the dotted lines are offset by factors 
of 2. 
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Fig. 8. Comparison of the XLF of S254-S258 (thick grey line) 
to the XLF of the O rion Nebula Cluster (from the COUP data; 
iGetman et al.ll2005bl black line). The straight lines show the re- 
sults of the power-law fits to the distributions in the luminosity 
range L x = 10 30 5 - 10 320 erg s -1 . 



(seelFeigelson & Ge tman 2001 IGetman et alJl2006HWang et all 



2007). 

To construct the XLF of S254-S258, we use the intrinsic full 
band [0.5 - 8.0] keV luminosities calculated by XPHOT (Table 3 
column (7)). Our XLF of S254-S258 is shown in Figure [8] and 
compared to the XLF for the stars in th e Orion Nebula Cluster 
from the COUP (IGetman et al.ll2005bl) . Obviously, the S254- 
S258 XLF peaks and turns over at a higher luminosity (near 
L x ~ 10 30 3 erg s _1 ) than the COUP XLF, because the X-ray 



detected sample for S254-S258 is incomplete for low-mass stars 
due to the lower sensitivity as discussed above. The slopes of the 
bright parts of these two distributions, however, can be seen to 
agree well. 

For a more quantitative analysis, we performed power-law 
fits of the form dN/d(log Lx) °c ax log L x for the observed dis- 
tribution of luminos ities. We used the maxim u m-lik elihood tech- 
nique described by Maschbe rger & Kroupal (l2009h . that yields 
an estimate for the exponent from the observed distribution func- 
tion (i.e. not a fit of the histogram). The resulting power-law ex- 
ponents for the distribution of X-ray luminosities in the range 
L x = [10 30 - 5 - 10 32 ] erg/s are a = -0.95 + 0.09 for the Orion 
COUP data, and a = -0.91 + 0.10 for our S254-S258 data. This 
consistency confirms t he results from comp arisons of other re- 
gions (e.g., the CCCP; lFeigelson et alj[201lh . 

This result shows that it is reasonable to assume that the 
XLF of S254-S258 has a similar shape as the ONC XLF. We 
can therefore make a quantitative estimate of the size of the to- 
tal young stellar population in the observed part of S254-S258 
by determining the vertical offset between the two distributions. 
We find that the total population in the observed part of the 
S254-S258 complex is m 0.7x of that in the ONC. Since the 
total population of the Orion Nebula Cluster (within 2.06 pc, 
Hillenbra nd & Hartmannl 19981) is about 2800 stars, the observed 
regiorfl of the S254-S258 should contain ~ 2000 stars in total. 

4.3. Spatial distribution of the X-ray sources 

The spatial distribution of the 364 X-ray sources in the Chandra 
field shows a complex pattern. Besides the prominent and dense 
central cluster S255-IR a few further apparent clusterings as 
well as a widely distributed population of X-ray sources that 
are spread homogeneously over the entire field-of-view of our 
Chandra observation can be seen. For a quantitative charac- 
terization of the spati al distribution we perfo rmed a nearest 
neighbor analysis (see ICasertano & Hutlll985l) to identify sta- 
tistically significant clusterings in an objective way. The surface 
density estimator at the position of each source i is given by 
AO'O) = _ 1)/ ( 7r -Dj(0 2 ) > where Dj(i) is the angular distance 
from source i to its jth nearest neighbor. We used 7 = 5 for our 
analysis; this value is large enough to keep the fluctuations of 
the local density estimates reasonably low and at the same time 
allows to detect groups with a minimum of ~ 5 stars. For the 
interpretation of the resulting densities, we have to take the spa- 
tial variations of the detection sensitivity over the field-of-view 
into account. The sensitivity is highest in the center, but due to 
effects such as mirror vignetting and the increasing width of the 
point-spread function, it gets several times worse near the edges 
of the ACIS field-of-view. We therefore plot in Fig. [9] the source 
density [x$ as a function of the offaxis-angle. The general trend 
of decreasing source density with increasing offaxis-angle can 
be clearly seen. 

Clusters can be defined as spatially confined groups of 
sources for which the local surface density clearly exceeds the 
values found at other locations in the image at similar offaxis- 
angles. The obvious dense central cluster S255-IR appears (as 
expected) as a very prominent peak at low offaxis-angles in 
Fig. [9] Using a density threshold of ^5 > 22 arcmin 1 , we find 
that 45 X-ray sources can be considered as members of this clus- 
ter. Two further prominent peaks can be seen in the plot: one 
peak consisting of 12 sources near offaxis-angle 5' is caused 



7 For comparison, we note that the diameter of the Orion Nebula 
Cluster (as 30") would be « 8' at the distance of 1.6 kpc. 
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Fig. 9. Top: Nearest neighbor analysis surface density at the lo- 
cation of each X-ray source plotted against the offaxis-angle. 
The general decrease of density with increasing offaxis-angle 
is related to instrumental effects. The dot-lined boxes mark the 
members of the three identified clusters. 

Bottom: Spatial distribution of the X-ray sources. The members 
of the three clusters marked in the top plot are show by the col- 
ored asterisks. 



5. The size of the distributed X-ray population 

In order to check whether the distributed X-ray sources may rep- 
resent the low-mass stellar population associated to the B0 stars 
in S255 and S257, we compare the size of the distributed X-ray 
population with the expe cted number o f based on IMF extrapo- 
lations. According to the iRroupal ( 200 1 ) field IMF model, each 
B0 star {Ma 15 M ; see iMartins et alj|2005T) should be associ- 
ated by ss 320 low-mass stars ([0.1 - 2] M ). However, not all of 
these low-mass stars will be detected in our X-ray observation, 
as the X-ray luminosities of young stars are related to the stellar 
mass. The expected number of X-ray sources can be found by 
comparing the X-ray detection limit to the typical X-ray lumi- 
nosity functions for young stars in specific mass ranges. 

As discussed above, the X-ray detection limit of our data 
is Lx.min ~ 10 29 5 erg s _1 in the central part of the observed 
field. However, this limit shows considerable systematic vari- 
ations as a function of the off-axis angle and gets several times 
worse near the edges of the ACIS field. If we consider the widely 
distributed population of X-ray sources, we have to take into 
account that most of these sources are located outsid e the cen- 
tral few arcmininute region of maximum sensitivity. iBroos et al.l 
(1201 lbl) performed a detailed analysis of the spatial sensitivity 
variations over the ACIS field. From the values for the com- 
pleteness limits in different off-axis angle slices in their Table 
8 we find that the area-weighted average of the completeness 
limit over the full field-of-view is about 0.5 dex higher than the 
on-axis value. This implies that the average X-ray completeness 
limit of our Chandra data for the widely distributed population 
is thus logLx ~ 10 30 erg s _1 . 

Since the X-ray luminosity functi ons for young stars are very 
similar for most stu died regions (see Feigelson & Getman 2005 



iGetman et al.l2006h we can assume that the young stars in S254 
S258 follow the same relations between stellar mass and X-ray 
luminosity as establishe d by the data from the Chan dra Orion 
Ultradeep Project (see iPreibisch & Feigelson! 12005). This im- 
plies that we should detect a 70% of the young stars in the mass 
range [0.5 - 2] M and « 30% of the stars in the mass range 

[0.1- 0.5] M Q 

Since the iKroupal (1200 ll) IMF predicts « 80 stars with 
[0.5 - 2] M Q and * 250 stars with [0.1 - 0.5] M Q for each B0 
star, the expected number of X-ray detected stars associated to 
the two B0 stars is ~ 2x(0.7x80+0.3x250) = 262. The observed 
number of as 250 distributed X-ray sources (after correction for 
background contamination) is actually quite close to this expec- 
tation value and consistent with the assumption that these X-ray 
sources trace the low-mass stellar population associated to the 
two B-stars. 



by a cluster of sources in the S256 region, while another peak 
around offaxis-angle 9' with 7 sources represents a clustering in 
the S258 region. 

Our clustering analysis thus reveals three significant clusters, 
that contain a total population of 64 X-ray sources. The remain- 
ing 300 X-ray sources are thus in a distributed, non-clustered 
spatial configuration. As discussed above, up to » 50 X-ray 
sources may be unrelated contaminants. This leads to a popu- 
lation of ~ 250 widely distributed X-ray detected young stars in 
the observed area. 



6. Conclusions 

Our deep Chandra observation of the S254-S258 complex led to 
the detection of 364 X-ray sources, about 50 of which are ex- 
pected to be background contaminants. The X-ray properties of 
most sources (luminosity, plasma temperature, and variability) 
are in the typical ranges found for young stellar objects. This 
supports the assumption that these X-ray sources objects repre- 
sent the population of young low-mass stars in the S254-S258 
complex. Our analysis of the spatial distribution of the X-ray 
sources with a nearest neighbor method reveals three signifi- 
cant clusters: the dense central cluster S255-IR, and two smaller 
clusterings related to S256 and S258. About 20% of the X- 
ray sources are members of one of these clusters, whereas the 
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large majority (~ 80%) of the X-ray sources traces a widely dis- 
tributed population of young stars. 

The size of this distributed population of X-ray detected 
young stars is in good agreement with the expected X-ray source 
number based on the assumption that these stars trace the low- 
mass population associated with the two early B-type stars in 
S255 and S257. We would not expect to see this distributed pop- 
ulation in the context of the models that two BO stars have either 
formed in isolation or were ejected from the central embedded 
cluster. Our results therefore suggest that the two B-stars and the 
associated distributed low-mass stars represent a stellar popula- 
tion that is distinct from the embedded cluster of YSOs in S255- 
IR. This is in agreement with the model scenario in which the 
observed star formation activity in the dense embedded cluster 
located in the interaction zone between the S255 and S257 H II 
regions has been triggered by the compression of the cloud due 
to the expansion of the H II regions. 

A detailed analysis of the optical and infrared properties of 
the individual X-ray detected young stars that can provide di- 
rect information on the ages, masses, and the circumstellar disks 
around these stars will be presented in a upcoming study. 

Acknowledgements. This work is based on observations obtained with 
the Chandra X-ray Observatory, which is operated by the Smithsonian 
Astrophysical Observatory for and on behalf of the National Aeronautics 
Space Administration (NASA) under contract NAS8-03060. Our analysis pre- 
sented in this paper was supported by the Munich Cluster of Excellence: 
"Origin and Structure of the Universe". This publication makes use of data 
products from the Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National Aeronautics 
and Space Administration and the National Science Foundation, and of obser- 
vations made with the Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology under a contract with 
NASA. 



References 

Allen, L. E., Hora, J. L., Megeath, S. T., et al. 2005, in IAU Symposium, Vol. 227, 

Massive Star Birth: A Crossroads of Astrophysics, ed. R. Cesaroni, M. Felli, 

E. Churchwell, & M. Walmsley, 352-357 
Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 
Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference 

Series, Vol. 101, Astronomical Data Analysis Software and Systems V ed. 

G. H. Jacoby & J. Barnes, 17-+ 
Bieging, J. H., Peters, W. L., Vila Vilaro, B., Schlottman, K., & Kulesa, C. 2009, 

AJ, 138, 975 

Brand, J., Massi, F., Zavagno, A., Deharveng, L., & Lefloch, B. 2011, A&A, 
527, A62+ 

Briceno, C, Preibisch, T., Sherry, W. H., et al. 2007, Protostars and Planets V 
345 

Broos, P. S., Feigelson, E. D., Townsley, L. K., et al. 2007, ApJS, 169, 353 
Broos, PS., Getman, K. V., Povich, M. S., et al. 201 la, ApJS, 194, 4 
Broos, P. S., Townsley, L. K., Feigelson, E. D., et al. 2010, ApJ, 714, 1582 
Broos, P. S., Townsley, L. K., Feigelson, E. D., et al. 201 lb, ApJS, 194, 2 
Casertano, S. & Hut, P. 1985, ApJ, 298, 80 
Cash, W. 1979, ApJ, 228, 939 

Chavarrfa, L. A., Allen, L. E., Hora, J. L., Brunt, C. M., & Fazio, G. G. 2008, 
ApJ, 682, 445 

Chopinet, M., Deharveng-Baudel, L., & Lortet-Zuckermann, M. C. 1974, A&A, 
30, 233 

Cyganowski, C. J., Brogan, C. L., & Hunter, T. R. 2007, AJ, 134, 346 
Deharveng, L., Zavagno, A., Schuller, E, et al. 2009, A&A, 496, 177 
Di Francesco, J., Johnstone, D., Kirk, H., MacKenzie, T., & Ledwosinska, E. 

2008, ApJS, 175, 277 
Favata, F, Flaccomio, E., Reale, F, et al. 2005, ApJS, 160, 469 
Feigelson, E., Townsley, L., Gudel, M., & Stassun, K. 2007, Protostars and 

Planets V, 313 

Feigelson, E. D. & Getman, K. V. 2005, ArXiv Astrophysics e-prints 
Feigelson, E. D., Getman, K. V, Townsley, L. K., et al. 2011, ApJS, 194, 9 
Forbrich, J. & Preibisch, T. 2007, A&A, 475, 959 

Freeman, P. E., Kashyap, V, Rosner, R., & Lamb, D. Q. 2002, ApJS, 138, 185 



Getman, K. V, Feigelson, E. D., Broos, P. S., Micela, G., & Garmire, G. P. 2008, 
ApJ, 688,418 

Getman, K. V., Feigelson, E. D., Broos, P. S., Townsley, L. K., & Garmire, G. P. 

2010, ApJ, 708, 1760 
Getman, K. V, Feigelson, E. D., Grosso, N., et al. 2005a, ApJS, 160, 353 
Getman, K. V, Feigelson, E. D., Townsley, L., et al. 2006, ApJS, 163, 306 
Getman, K. V, Flaccomio, E., Broos, P. S., et al. 2005b, ApJS, 160, 319 
Goddi, C, Moscadelli, L., Sanna, A., Cesaroni, R., & Minier, V. 2007, A&A, 

461, 1027 

Gudel, M., Briggs, K. R., Arzner, K., et al. 2007, A&A, 468, 353 

Heyer, M. H., Snell, R. L., Morgan, J., & Schloerb, F. P. 1989, ApJ, 346, 220 

Hillenbrand, L. A. & Hartmann, L. W. 1998, ApJ, 492, 540 

Howard, E. M., Pipher, J. L., & Forrest, W. J. 1997, ApJ, 481, 327 

Imanishi, K., Koyama, K., & Tsuboi, Y. 2001, ApJ, 557, 747 

Itoh, Y, Tamura, M., Suto, H., et al. 2001, PASJ, 53, 495 

Jaffe, D. T, Davidson, J. A., Dragovan, M., & Hildebrand, R. H. 1984, ApJ, 284, 
637 

Jiang, Z., Tamura, M., Hoare, M. G., et al. 2008, ApJ, 673, L175 
Kroupa, P. 2001, MNRAS, 322, 231 

Kurtz, S., Churchwell, E., & Wood, D. O. S. 1994, ApJS, 91, 659 
Kurtz, S., Hofner, P., & Alvarez, C. V. 2004, ApJS, 155, 149 
Longmore, S. N, Burton, M. G., Minier, V., & Walsh, A. J. 2006, MNRAS, 369, 
1196 

Martins, F, Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 

Maschberger, T. & Kroupa, P. 2009, MNRAS, 395, 931 

Minier, V., Burton, M. G., Hill, T, et al. 2005, A&A, 429, 945 

Minier, V., Peretto, N., Longmore, S. N., et al. 2007, in IAU Symposium, Vol. 

237, IAU Symposium, ed. B. G. Elmegreen & J. Palous, 160-164 
Miralles, M. P., Salas, L., Cruz-Gonzalez, I., & Kurtz, S. 1997, ApJ, 488, 749 
Moffat, A. F. J., Jackson, P. D., & Fitzgerald, M. P. 1979, A&AS, 38, 197 
Ojha, D., Tamura, M., & Sirius Team. 2006, Bulletin of the Astronomical Society 

of India, 34, 119 
Pismis, P. & Hasse, I. 1976, Ap&SS, 45, 79 
Preibisch, T. & Feigelson, E. D. 2005, ApJS, 160, 390 
Preibisch, T., Hodgkin, S., Irwin, M., et al. 2011, ApJS, 194, 10 
Preibisch, T., Kim, Y, Favata, E, et al. 2005, ApJS, 160, 401 
Preibisch, T. & Zinnecker, H. 2002, AJ, 123, 1613 

Preibisch, T., Zinnecker, H., & Schmitt, J. H. M. M. 1993, A&A, 279, L33 
Rygl, K. L. J., Brunthaler, A., Reid, M. J., et al. 2010, A&A, 511, A2+ 
Sharpless, S. 1959, ApJS, 4, 257 

Smith, R. K, Brickhouse, N. S., Liedahl, D. A., & Raymond, J. C. 2001, ApJ, 
556, L91 

Snell, R. L. & Bally, J. 1986, ApJ, 303, 683 

Stassun, K. G., van den Berg, M., Feigelson, E., & Flaccomio, E. 2006, ApJ, 
649, 914 

Testi, L., Palla, F, & Natta, A. 1999, A&A, 342, 515 
Townsley, L. K, Broos, P. S., Corcoran, M. F, et al. 2011, ApJS, 194, 1 
Townsley, L. K, Feigelson, E. D., Montmerle, T., et al. 2003, ApJ, 593, 874 
Wachter, K, Leach, R., & Kellogg, E. 1979, ApJ, 230, 274 
Wang, J., Townsley, L. K, Feigelson, E. D., et al. 2007, ApJS, 168, 100 
Wang, Y, Beuther, H., Bik, A., et al. 2011, A&A, 527, A32+ 
Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 
Wolk, S. J., Harnden, Jr., F. R., Flaccomio, E., et al. 2005, ApJS, 160, 423 
Zavagno, A., Russeil, D., Motte, F, et al. 2010, A&A, 518, L81 + 
Zinnecker, H., McCaughrean, M. J., & Wilking, B. A. 1993, in Protostars and 
Planets III, ed. E. H. Levy & J. I. Lunine, 429^195 



